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Polarized entangled Bose-Einstein condensation 
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We consider a mixture of two distinct species of atoms of pseudospin-^ with both intraspecies and 
Interspecies spin-exchange interactions, and find all the ground stats in a general case of the param- 
eters in the effective Hamiltonian. In general, corresponding to the two species and two pseudo-spin 
states, there are four orbital wave functions into which the atoms condense. We find that in certain 
parameter regimes, the ground state is the so-called polarized entangled Bose-Einstein condensation, 
i.e. in addition to condensation of interspecies singlet pairs, there are unpaired atoms with spins 
polarized in the same direction. The interspecies entanglement and polarization significantly affect 
the generalized Gross-Pitaevskii equations governing the four orbital wave functions into which the 
atoms condense, as an interesting interplay between spin and orbital degrees of freedom. 

PACS numbers: 03.75.Mn, 03.75.Gg 



I. INTRODUCTION 

Since a decade ago, there have been a lot of activi- 
ties on Bose-Einstein condensation (BEC) of atoms with 

fin degree of freedom [l[, for example, in spin-1 0- 
and pseudospin-i gases [7l-[l0j. as well as BEC of a 
mixture of two kinds of spinless atoms [U, [ll| - [T7| . As a 
step beyond these research hnes, a novel class of BEC, 
the so-called entangled BEC (EBEC), was proposed as 
the ground state of a mixture of two distinct species of 
pseudospin-i atoms with interspecies spin exchange in 



some parameter regime |18l-l20|. Spin-exchange scatter- 
ing between atoms of different species entangles these two 
species, and EBEC may lead to a novel kind of super- 

I 



fluidity [2l|. More recently, EBEC has also been found 
in a mixture of spin-1 atoms [22|. It has been noted 
that EBEC may be experimentally realized using two 
different species of alkali atoms. In some aspects, EBEC 
bears some analogies with a single species of pseudospin- 
i atoms in a double well or occupying two orbital 
modes [1, 0] , but there are also important differences due 
to the fact that two atoms of distinct species are distin- 
guishable. In general, there are four orbital wave func- 
tions into which the atoms of the two species and two 
pseudospin states condense into, respectively. 

Under single-mode approximation for orbital degree 
of freedom, the many-body Hamiltonian can be written 
as 



(1) 



where aa (a = a,b, a 
of the species a, 



=t, i) is the annihilation operator 



I 

^-Kh? ^'"^^)^, ^ / ^ap, where a and a' may or may not be 
equal. For interspecies scattering with spin exchange, 



X 
J3 



K 



(a/3) 



is related to the scattering lengths as the following. In 
general, ^j^^lasa^ is the scattering length for the scat- 
tering in which an a-atom flips from cr4 to ai while a 
/3-atom flips from 0-3 to (T2, fiai3 = mamp/^ma + rrifj) 
is the reduced mass. Now, for intraspecies scattering, 
(7ctct"'' = 27r?i^^CTCTCTcr/Mcta for same species and same spin. 



For interspecies scattering without spin exchange, 



(ab) 



5e / 0L(r)^^^(r)0b,(r)</),^(r)d3r, (3) 



where a ^ a,ge = giflg = 2-!:%^ S}^"^^ / ^i^p ■ faa = ^aa- 
K'ff^ jl, where t^a = / (i>*aa^aa'\>aad^-r is the single- 
particle energy, h^ry being the single-particle Hamilto- 
nian. 

The total spin operator of species a (a — a, b) is 
Sq ~ a^^Sa-a'Ct'g., where Sa-a' is the single spin operator. 
Hence the Hamiltonian can be transformed into that of 
two coupled giant spins. 
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— '2'Ke[SaxSbx + SaySby) + JzSazSbz + BaSaz + B^Sbz + ^aS'^z 



Eq, 



(4) 



where Jz 



X{f ^-Xjf \ Ba 



N^(T^{a.a.) T^(o-o-)\ , N,, ( T^(a.b) 



K. 



(ab) 



2 v-"-n ^ 
= At - hi + 1 - Kf^^) + ^{K. 



/at-/f4+ 

I ^ I I ^7 



Aab) 
(aa) 



(bb) 



.{bb). 



it 



2 V'^'^tt 

(l/2)Ea(E./o.)^a + (1/8)E.(E..'4:'V^ + 

(1/4) (Eaa' Klf})NaNb doBS xiot depend on spins. We 
label the two species in such a way that Na > Nb- 

As discussed previously, at the antifcrroniagnctic 
isotropic point of the effective parameters, namely 2Ke = 
Jz > while Ba = Bb = Ca = Cb = 0, n = 
2KeSa ■ Sb. Then the ground states are \ ^\^\ Sz) = 
^(„t)JV„/2-JV,/2+s.(^piv„/2-iv./2-s.(^t;,t _ a^fep^'lO), 
where Z is the normalization constant. This ground state 
leads to interesting physical consequences. It has been 
shown numerically and analytically that EBEC also per- 
sists in a considerable regime away from the isotropic 
point. 

In this paper, we extend the consideration of this 
model to a more general case of the parameters, namely, 
Ba = Bb, denoted as B, Ca = Cb, denoted as C, and 
Jz — 2Ke — 2C — 0, and thus the Hamiltonian becomes 



H = 2KeS^ ■ Sb + BSz + est + Eq. 



(5) 



I 

the form of \S, Sz)- The ground state is thus 

|G>-|5",5r), (6) 

where 5™ and 5™ are, respectively, the values of 5 and 
5, that minimize 



E ^ KeSiS + 1) + BSz + CSl 
In case C 7^ 0, it can be rewritten as 



(7) 



E = K,(^S+^^ +c{Sz + ^^ +Ei, (8) 

where Eq is independent of S and Sz ■ The minimization 
should be under the constraints 



Q . <<?<<? 

>~^rmn _ ^ _ ^maxi 

-S <Sz< s, 



where 



Smax — Sa ~l~ Sh — 



Na~Nb 



Na + Nb 



(9) 
(10) 



(11) 



(12) 



which reduces to the isotropic antiferromagnetic coupling 
when B = C — while Ke > 0. In other words, we 
assume that there is such a rotational symmetry that 
the Hamiltonian only depends on the total S and Sz, 
which are thus good quantum numbers. We believe that 
under a perturbation due to deviation of the parameters 
from this condition, the ground state is still close to the 
present one in each parameter regime discussed below, 
based on an extension of the previous analysis for the 
isotropic parameter point j20| . the details of which will 
be reported elsewhere. 

In Secini we find all the ground states of the Hamilto- 
nian ([SJ in various parameter regimes, based on the cal- 
culation accounted in the Appendices. In Sec. IIIIl these 
ground states are expressed in terms of Bosonic degrees 
of freedom. Most of the ground states are of the form 
of 15*, its'), which are referred to as the polarized BEG 
if S' ^ 0, as then 5*2 = ±5* ^ 0. Their properties are 
discussed in Sec |lVl A summary is given in Sec. El 



II. GROUND STATES IN VARIOUS 
PARAMETER REGIMES 

For spin-i atoms, Sa = Na/2 and Sb = Nb/2 are fixed. 
It is clear that the total spin S and its z-component Sz 
are good quantum numbers, hence the eigenstates are of 



We have obtained all the ground states in the whole 
three-dimensional C — Ke — B parameter space. 



A. B / 

Based on the calculations in Appendix lAl the ground 
states for i? 7^ are summarized in Table U and in the 
two-dimensional C — Kf. phase diagrams for a given B 7^ 
0, as shown in FIG. [TJ as well as FIG. [5] for the special 
case of Na^ Nb^ N. 

In the regime < ^\+f for C < plus < 



\B\-{Na+N^)C 



N^+N,+i for < C < ivlfk' ground state 

is |^^^^,-sgn(B)(^^^^±^)), where sgn(B) denotes the 
sign of B. 

In the regime > ^\'^f for C < plus > 

^''^NMi'' for < C < ^J^, the ground state 

is |^^^^,-sgn(B)(^^^^)). Note that if Na - Nb, 

I B I 

then this statement is still valid by regarding as 



infinity. In other words, the boundary C 



\B\ 



disap- 



pears, and the ground state is |0, 0) in the whole regime of 
' \+f fo^ <-\B\ plus Ke > \B\ioiO-\B\. 



Ke > 



On the boundary Ke — ^^L ^"'^ for 



C 



< 



\B\ 



^i^±^,-agn{B){^H^)) and 



'-)) 
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Parameter regimes 


Ground states |S"",Sr> 


B / 


C<-\B\ 






sgn(B)(^"r")) 


_ IBI-NaC 




^,_sgn(B)(M))and 
sgn(S)(^»r^)> 


^ N„.+l 






C = -\B\ 


Ke < -C 




--r^, sgn(B)(-»r^)> 


K, = -C 




5"", sgn(B)5'"), ^--^>- KS"" < '^"l^'' 


> -C 




^^,-sgn(i3)(^^)> 




1 D 1 / 7V7" 1 AT" ^ 




AT" 1 AT" AT" I AT 

^^^^,-sgn(B)(^^^^^)) 


1 D 1 I AT 1 AT \ 1 D 1 / AT AT \ 

\D\-(No,-^Nk)0 ^ r<' ^ l-Hl-CJVa-JVbJC 










s, sgn(B)s), s - Int(J(^^_^^)) 


IV ^ \B\-{Na+Nf,)C 

- iVa+JVh + 1 




iY^,_sgn(B)(iX^)> 


1 R 1 1 Rl 


-fSTe <0 






/^e =0 




5'",Int(-^)), llnt(-^)| < 5"^ < 


1 Rl ^ jv j\r, 




... / Lb 1 Fc \ 


Na+N^ - - JV„-JV(, 




s, sgn(B)s), s - Int(2'(^^_^(5,) 


tr \B\-(N^--N^)C 
" — N^-Nh + l 




sgn(S)---^) 


r ^ \B\ 


-fSTe <0 




™,Int(-^)> 


/STe =0 




5-,Int(-^)), |Int(-^)| < 5- < 






> 




^,Int(-#.)> 



TABLE I. Ground states in various parameter regimes with B ^ 0. 



{ '^"2'^'' , -sgn(-B) ( ^° 2 )) ^"^^ degenerate 
ground states. 

In the regime < k < W-(^<'-^^)c 

and ife > for -\B\ < C < jj^^, the ground 
state is |s, — sgn(i?)s), -where s = Int(S'o) is the inte- 
ger closest to So = 2{K^+c) ^^'^ legitimate range 
Smin < s < Smax- When Na = Nb, the boundary 

^'^ = '^'v^j^iVT+i^^ becomes K,, = \B\ for C > -\B\ 
up to infinity. Tliis is a crossover regime, as s depends 
on the parameters, and is continuously connected with 
the two neighboring regimes, that is, s = at the 



boundary = ^^If+Jl^"" while s = 
boundary K, 

ities on the boundary Kg = for 



• "^^^ there are discontinu- 



In the regime C > 



I-BI 



state is 



\Na+Nb 



while ii'e < 0, the ground 



-Int(2^)). It is continuously connected 



with |iXa±iVL^_sgn(B)(^^4^)) on the boundary C 



\B\ 



while <0. 



-sgn(B)(^)) 



Na+Nb 

I S I 

Finally, for Na > Nh, in the regime C > while 
Ke > 0, the ground state is 1^^^=^^, Int(-^)). It is 
continuously connected with | , 

I B I 

on the boundary C = while Kf, > 0. 

I-BI I-BI 

On the boundary = while < C < , 

the degenerate ground states are Int(— ^)), with 
|Int(-^)| < 5™ < On the boundary = 



2C 

while C > 



\B\ 



N„-Nb ' 



the degenerate ground states are 
|5'",Int(-^)), with ^i^^ < S"" < The 
ground state in each regime bordering each of these two 



boundaries is one of the ground states on it. Hence dis- 
continuities or quantum phase transitions occur in cross- 
ing this boundary from one regime to another. 

The only other boundary where discontim.iity or quan- 



tum phase transitions occur is Ke 



\B\-N„C 

Na + 1 



for C < 



— \B\. In entering it from downside or upside, the 
ground state remains as the original one, but when leav- 
ing it and entering the other side, the ground state 
discontinuously changes to another state. Neverthe- 
less, if entering the boundary from the crossover regime 

jy^f , there is no discontinuity, as 



-\B\ < C < 



the groimd state \s, —sgn(B)s) in this regime is contin- 



I Ng+Nb 

2 



-sgn(S)( 



)) 



uously connected with both 
and|M,_sgn(i3)(M)). 

Three boundaries Kg = ^^j^^l^ for C < —\B\, = 
and = ^^^^ for C > -\B\ con- 
verge at Ke = —C = \B\, which does not belong to any 
of these three regimes and where the ground state can be 
any state \S"',S^), with S"" and in the legitimate 
ranges, or any superposition of these states. In the three 
dimensional C — — B parameter space, these three 
boundaries that arc lines on the C — plane with a 
given B become planes. = —C = \B\ are two lines 
converging at the origin for B > and B < Q respec- 
tively. 

The structure of the three-dimensional phase diagram 
possesses a reflection symmetry with respect to the plane 
B = Q. In each pair of symmetric regimes, the ground 
states have the same S*™ but opposite 5™. 

As iJ — >■ 0, the boundary C — for <{) ap- 



4 




FIG. 1. Ground states in C — Ke parameter plane for a given B ^ 0, which can be positive or negative. The solid lines are 
boundaries between different regimes. The dashed lines are only coordinate axes. The ground states in five regimes are indicated. 
|s, — sgn(i3)s), with s = Int( 2(^~+c) )i is the ground state in the regime surrounded by Ke = "^^j^^'j^'^^j'"''^ , Ke — 
and Ke = , which is continuously connected with | '^"'^^i' ^ — sgn(i3)( )) and | '^'^~^>' ^ —sgn{B){ ^'^^^'' )) on the two 

boundaries, respectively. The boundary for C < —\B\ is Ke = ^—j^-^^, on which both | , — sgn(_B)( )) and 

I , —sgn{B){ '^°~^'' )) are degenerate ground states. Three boundary lines converge at Ke = —C — \B\, on which any 

state is the ground state. The ground state '^"■^^f> ^ _lnt(^)) is continuously connected with \!^-^^^, — sgn(_B)(-^^^-i^)) on 
the boundary C = ^ ^"^^^ while Ke < 0. Int( — ^)) is continuously connected with | '^'^~'^>' ^ ~sgn{B){ ^'^~^'' )) on the 

\ B\ \ B\ \ B[ 

boundary C — -m^ while Ke > 0. On the boundary Ke ~ while ' < C < ^ , the degenerate ground states are 
of the form of |S"\Int(-^)) with |Int(-^)| < S"" < '^''+^» . On the boundary Ke ^ while C > jf^^, the degenerate 
ground states are of the form of |S"", Int(— 2^)) with < 5"" < '^'^+^b ^ -pj^jg figure reduces to FIG.[2]in the special case 

0iNa=Nt. 



proaches the negative half of Kg-axis including the ori- 
gin, while the boundary C = ^ for Kg > 0, exist- 
ing when Na > Nt, approaches the positive half of Ke- 
axis excluding the origin. The boundary Ke — ^^l^^l^ 
for C < -\B\ approaches Ke = for C < 

, while the crossover regime of |s, — sgn(i?)s) tends to 
vanish, and Ke = — C = \B\ approaches the origin 
Ke = -C = \B\ = 0. 

B. 5 = 

Based on the calculations in Appendix |Bl the ground 
states for i? = are summarized in Table HIl and depicted 
in the two-dimensional C — Ke phase diagrams for for 
B = 0, as shown in FIG. d 

For C > 0, i^e > while B 0, the ground state 



is uniquely \—^ — ^,0). In the case of Na > Nh, it is 
continuously connected with | ^ int(— ^)) in the 

I B I 

regime C > _ , Ke > while i? 7^ 0. In the case 
of Na — Nb, the ground state for C > 0, Ke > while 
B — becomes |0, 0), which is the same one as for C > 
while Ke > \B\. 

For C > 0, Ke < while B — 0, the ground state is 
uniquely | ^''+^'' ^ 0). It is continuously connected with 

\^i^,-lnti^)) in the regime C > Ke < 

while B ^0. 

On the boundary C > 0, Ke = while B = 0, i.e. 
on the positive C axis, the degenerate ground states are 
|S"",0), with ^^^^ < 5™ < ^^4^, which include the 
two ground states in the regimes above and below this 
half axis. In case Na > Nf,, these degenerate ground 
states on positive C axis is continuously connected with 
|5'",Int(-^)), with ^^^^ < 5" < on the 



5 




FIG. 2. Ground states in C — Ke. parameter plane for a given B 7^ in the special case of A'^a = Nt — N . This diagram is a 
limit of FIG. \T\ with the boundary C — ^"^^^ disappearing while the boundary Ke = ^^^^^^"j^^^^^^^ becoming a horizontal 
line A'e = \B\. 



Parameter regimes 


Ground states \S"\ ST) 


B = 


C > 


Ke>0 




Ke = 




Ke<0 




C = 


Ke>0 




Ke = 




Ke<0 


1 :j 1 -Jz ^ "Jz ^ ^ 


C < 


> ]V„+1 




ry- _ iVaC 
-^■^ - N„+l 


l^^-^Si^"-^") and r''+^^±'^»^^^) 


< ]V„+1 


2 ' ^ 2 / 



TABLE II. Ground states in various parameter regimes with B = 0. 



boundary = 0, C > while B ^ 0. In case 

Na = A^f), these ground states are continuously connected 
with the ground states \S'"\ Int(-^)), with Int(-^) < 
g,n < E,l±Ej>^ on the boundary ^ 0, C > 
while B ^0. 

For C < 0, Ke > while B = 0, the twofold 

degenerate ground states are | ^ ± ). Simi- 

larly, for C < 0, ife < while B = 0, the twofold 

degenerate ground states are | ^°+^'' , zfc ^°+^'' ). On the 
boundary C < 0, ifg = while 5 = 0, the four- 

fold degenerate ground states are | ^ ± ^a-N^ ^^ 



I ^°+^'' , zfc ^''+^'' ), just the union of the ground states 
above and below this axis. 

For B = 0, on axis, the ground states are always 
degenerate too. For C — 0, Kf. > while B = 0, i.e. on 
the positive Kg axis, the ground states are | , 5'™), 

with -i^i^ < S-™ < i^^^^, which include the ground 
states in the regimes on the right and left of this half 
axis. This is the ground states discussed previously. For 
B — 0, C — while Ke < 0, i.e. on the negative K^ 
axis, the ground states are S^}, with - ^°+^'' < 

grti £ Na+Nb ^ ^]^[(.]^^ again, include the ground states in 
the regimes on the right and left of this half axis. 
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0) 



c 



2 ' 



± 



,0) 



FIG. 3. Ground states in C — Ke parameter plane for B — 0. The solid lines are boundaries between different regimes. 
The dashed line is only the negative C axis. The ground states in four regimes are indicated. The boundary for C < is 
= on which both | ^°+^% ± ^°+^'' ) and | ^ -b ^--^" ) are degenerate ground states. On the positive C axis, 

the ground states are |S"\0>, with < S"" < ^°+^\ On the positive Ke axis, the ground states are l ^"'^" , ST) , 

with < ST < ^"2^" ■ On the negative axis, the ground states are \^^^^±^,ST), with < 5™ < iX^L^Xk. 

On the origin B = C ^ = 0, the ground state can be any superposition of 15"", SD, with < 5"" < and 



At the origin B = C = — 0, the ground state can 
be any legitimate state, i.e. |S"",S'^), with < 
5" < ^°+^'' and -S"" < S"^ < 5"", or any of their 
superposition. 

Therefore in the case of i? = 0, all the boundaries on 
C — Ke parameter plane is a phase boundary, on which 
discontinuities or quantum phase transitions occur. The 
degenerate ground states on the boundary include the 
ones in the regimes it divides. 

As i? — > 0, the boundaries of different C — Kf, regime 
for 7^ converge to those of B = 0. A ground state 
in each C — Ke regime on each side of S = reduces to 
a ground state in the corresponding C — K^ regime for 
B = Q. For C > 0, all the ground states for i? 7^ are 
continuously connected with those for B = Q. For B = {) 
while Ke < 0, the ground states in the regime C > can 
be obtained by substituting i? = in the ground state 
for C > jvlfivr while Ke < 0. For B = while Ke > 0, 
the ground states in the regime C > can be obtained by 

I B I 

substituting i? = in the ground state for C > 

while Ke > in case Na > Ni,, and it is the same as the 

ground state for C > while Ke > \B\ in case Na = Ni,. 



For C < while B = 0, excluding the origin B = 
C — Ke = 0, the ground states in each C — Ke regime 
include ground states on both sides of i? = 0. Hence for 
C < 0, there is always a discontinuity or quantum phase 
transition in the ground state when B switches its sign 
from positive to negative and vice versa, for any values 
of Ke- When B continuously reduce the magnitude, the 
ground state remains as the initial one when reaching 
B = 0, then the discontinuity in ground state occurs 
when the sign B becomes opposite to the original. 

III. GROUND STATES IN TERMS OF BOSON 
OPERATORS 

One can obtain the expression of any ground state 
IS", Sz) in terms of bosonic degrees of freedom, using 




(13) 

where g{S,Sz,Shz) = {Sg.S^ - Sbz, Sb, Sbz\S, -S) is the 
Clebsch-Gordan coefficient 1231 
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{2S + l){Sa + Sb- S)l{Sa -Sb + S)l{Sb -Sa+ S)l 



9(^7 ^z, Shz) — {Sa, Sz Shz', Sb, Sbz\Sa, Sh, S, Sz) ~ +S+1)! 

X [{Sa + Sz- Sbz)l{Sa -Sz + Sbz)l{Sb + SbzV.{Sb ^ Sbz)l{S + - 5,)!]^ 

X Efe(-l)'[fc!(5a + Sb-S- k)l{Sa -Sz + Sbz - k)\{Sb + Sbz - -Sb + Sz- Sbz + k)\{S - Sa ~ Sbz + fc)!]-\ 

(14) 

where k is an integer such that the arguments in the factorials are non-negative. In terms of the particle number in 
each single particle mode a = a, 6, 



\Son Saz) a — | "5'q ~\~ Saz)a\\Sa Saz)a].-: 



(15) 



where \n)aa 
Therefore, 



^(at)"|0), {a =t,;). 



Si 



\s,Sz) - /(^,^.,^6.)(4)^''+^^'^-(4)^"-^^+^-(4)^^+^-(&l)^^-^-io), 



(16) 



Shz—~Sh 



with 



f{S, Sz, Sbz) = 9{S, Sz, Sbz)[{Sa +Sz- SbzV-iSa -Sz+ Sbz)KSb + Sbz)KSb - SbzV-]-^^^. 
In some special cases, to which most of our ground states belong, this expression can be simplified as the following. 



For Sz = —S, we have 

g{S,-S,Sbz)^{-lf'+'''^ 



Therefore 



where Fi = | 



fiS,-S,Sbz)^i-l) 



A. \S,±S) 

{2s + 1)!(5„ + Sb- sy.jSq + Sbz + sy.jSb - StzV. 
iSa + Sb + s + ly.iSa -Sb + sy.{Sa ~s- SbzyiSb + Sbzy 

{Sa + Sb- sy 



Sb-\-Sbz 



ri 



(Sa — s — Sbzyi^b + Sbzy ' 



(25+1)! 



1/2 



(S,„„x + S'+l)!(S„ax-S)!(S„i„ + S)! 



Thus 



\S,-S) = r,{alf+^{blf ^{a\bl-alb\)^'^\0). 



Similarly, for Sz = S, we have 



(25 + 1)!(5, + Sb- Sy.{Sa - Sbz + sy.{Sb + Sbzy 



Therefore 



Thus 



/(5,5,56,) = (-l) 



iSa + Sb + S+ l)!(5a -Sb + 5)!(5a - 5 + SbzyiSb - Sbzy 
{Sa + Sb-Sy 



Fi 



{Sa - s - SbzyiSb + Sbzy 



Wg + JVi, 



\S,S) ^Tiia\f+^^ib\f~^^ia\bl-alb\)^^^-''\0). 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



It is clear that in |5, ±5), there are — 5 interspecies singlet pairs, plus 5 + a-atoms and 5 — 

6-atoms, each of which has z-component spin polarized as ±1. This is fully consistent with discussions in terms of 
the generating function method |22| . 
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B. |^^,5.> 



For 5 — Smin , we have 



l(Smim Sz, Sf)z) — ( 1) 



i2Sa ~ 2Sb + iy.{2Sby.{Sa +Sz~ StzY-iSa -Sz+ SbzV. 



{2Sa + + SbzV-iSb - SbzV-iSa -Sb + SzY.iSa -Sb- SzV. 



Therefore 



f{Smim Sz, Sbz) — r2(— 1) 



iSb + Sbz)KSb-Sbzy- 



{(2S —2S +1)! 1 "^^^ 

(23^ + 1^"'^ ^ ^'^^ O M/nC M !• . Thus 



)!(S„-5b+S,)!(S„-Sb-S,)!(2St)! } 

,Sz)=T,{a\)^+'^ial) 



■■ia\bl-alb\Y 



as obtained in \li 



2 ' 



(23) 



(24) 



(25) 



For 5 = Smax, we have 

di^max, Sz, Sbz) 



{2Sb)K2Sa)KSa +Sb + SzY-iSa + Sb - S z)l 



Therefore 



\Ei±I^^^Sz) = r, c(5fc.)(4)^"+^^"^^^(4)^"-^^+^^^(4)^''+^''^(4)^''-^^10), 



{2Sa + 2Sb)KSa + Sz- SbzV-iSa -Sz + Sbz)KSb + Sbz)KSb - Sbz)l 



1/2 



Sbz——Sb 



where T, = { (^g-.)K2S,)!(g.+S^+S.)!(5,+g^ ^ c(5,,) = [{3^ + Sz ^ SbzYiSa - Sz + Sbz)l{Sb + SbzYiSb 

SbzY-]"^, which is nonzero only if the all the arguments of factorials are non- negative. 



(26) 



(27) 



Especially, when Sz = ^{Sa + Sb), the state is a direct product of states of individual species. 

Na±Nb ^ Ng + Nb ^ _ Na Na. .1% Nb. 

I 2 ' 2 ^ ~ ' 2 ' 2 2 ' 2 



where 



' 2 ' 2 



(a 



|0), 



(28) 



(29) 



with Oj. representing al and at, is a ferromagnetic state of species a. 



r 



IV. ENTANGLEMENT 



sider the entanglement between the occupation number 
of each of the four single-particle modes and the rest of 
the system Alternatively one simply considers the 
In any state \S, Sz), the two species are entangled un- total the spins of the two species. Both approaches yield 
less S = Sa + Sb and 5"^ = ±5. In general, we can the same entanglement entropy £, 
call a state polarized entangled BEC when Sz and 
there is interspecies entanglement. For any state \S, Sz), 
we calculate the entanglement between the two species, Sb 

which is quantified as the von Neumann entropy of the ^ = ^ X] 3^^^' ^°&Nb+i 9'^{Sz,m). (30) 

reduced density matrix of either species. One may con- m=-Sb 
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The reason why both approaches lead to the same re- 
sult, for any jS', S'^), is that the expansion is always 
automatically Schmidt decomposition with coefficients 
g{S, Sz,Tn), no matter it is in terms of the spin basis 
states of the two species, as given in (IT3)) . or in terms of 
the occupation numbers of each species and each pseu- 
dospin state, as obtained by substituting to (|13l) . 

In the following, we focus on the maximally polarized 
state IS", its'), because of the following reasons. First, 
most of our ground states are in this form. Second, in 
case the ground states are degenerate with respect to 
Sz , an infinitesimal symmetry breaking perturbation pro- 
portional to Sz picks out a maximally polarized state 
as the new ground state. Clearly, the entanglement 
entropies of \S,S) and of jS*, — 5) are equal, because 
g{S,m) = [-lf^"^-Sg{-S,-m). 

For simplicity, we reduce our focus to the case of 
Na ^ Ni, ^ N, in which \S,-S} = Ti{alblf {a\bl - 
46|)^^^|0) and \S,S) ^Ti{a\b\)^{a\bl~ alb\)^-^\0). 

The calculation results are shown in FIG. |4] for four 
values of N. It can be seen that for a given N, entangle- 
ment always decreases when S increases. 

We know 

iV„, = iV„/2 + ry,5a„ (31) 

with ?7-|- = 1 and r/^ = —1. Therefore, for any state, the 
I 



average number of a-particles with spin a is 

N 

{Nc.„)^-+f^,M^u (32) 

where 

M,i EE {Saz), (33) 

Its fluctuation is 

A7V„. ^ ^{Nl,) - ((iV„,))2 = (M„2 - Ml^fl\ (34) 
where 

M„2^(5L). (35) 

For \S,-qsS), with 77s = ±1, it is obtained that Mai — 
Mbi — r]s§. Hence the total polarization is rjsS. It is 
also obtained that for this state, Ma2 = Mb2, denoted as 
M2, which depends on S and A'^ in a more complicated 
manner, as shown in FIG. El For this state, one also 
has ANa^ — ANai — ANt^ = AiV^^, which is shown in 
FIG [51 One can see that the dependence of the fluctua- 
tion ANaa- on S has a trend similar to the entanglement 
entropy. Indeed, the fluctuation ATV^^ is also a charac- 
terization of the interspecies entanglement [19 ] . 

Both Mai and M2 enter the Gross-Pitaevskii-like equa- 
tion governing the single-particle orbital wave function 
4>aa-ir) under the ground state IS^ijsS), with S < 





5\ r 







2mn 



4 



{N~1)S 



M2-^ 



(36) 



It is obtained by requiring the variation of the energy 
with respect to (paa to be minimum, under the constraint 
/ \4'a<y\'^dv ~ 1. Note that now the energy {H), where T-L 
is as given in ([T]), is treated as a functional of (jiaa- The 
effective coefficients tea and if's are now functionals of 
4>aa- In limit S* = 0, the is equation reduces to that 
under the singlet ground state previously studied. 

It can be seen that interspecies entanglement sig- 
nificantly affects the Gross-Pitaevskii equations, noting 



I 

M2 = {ANaa? + ^V4. 

When S — N, the ground states \N,riNN) is disen- 
tangled between the two species, consequently there are 
only (paa- (oi = a,b) with a specified by rj^rjj^ = 1. For 

5| |0), there are only 0at and (j)h\- 



\N,N) ^ 
For \N,-N) 



(l/iV!)a| 



6| |0), there are only (pai and 
(ph^. For these two states, AI2 = iV^/4, and thus 
reduces to 



2m, 



+ (TV - l)g 



{aa)\ 



(37) 



which is just Gross-Pitaevskii equations for the well stud- 



ied two-component BEG. 
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FIG. 4. Entanglement entropy of the interspecies entanglement in state IS, izS), with A'^a — Nb = N , as a. function of S. 



V. SUMMARY 

To summarize, we have considered a more general case 
of the many-body Hamiltonian of a mixture of two species 
of pseudospin-i bosons with interspecies spin-exchange 
interaction, with three effective parameters as given in 
([5]). We have determined the ground states in all the 
regimes of these three parameters. Discontinuities of 
the ground states occur in some parameter boundaries. 
These are first order quantum phase transitions, as S"^ or 
have discontinuities in crossing the boundaries. 

The two species are entangled in all the ground states 
except | ^°^^'' ,ib ^''^^'' ), which is the ground states 
when K,, and C are both negative or sufficiently small 
positive, as shown in the phase diagrams. Very interest- 
ingly, when Na = Nb, the ground state is the global sin- 
glet state |0, 0) in the wide regime in which Ke is greater 
than the larger one between and \B\. This result 

confirms the previous arguments based on numerical esti- 
mation and perturbative analysis that the singlet ground 
state and the interspecies entanglement persist in a wide 
parameter regime [l^, [l^l . 

We have calculated the properties of \S, ±S) in details. 
For a given N, the polarization is equal to ztS. The 



entanglement decreases with as S increases. It reaches 
the maximal value 1 as iS = 0, and reaches the mini- 
mal value as 5 = A^. The fiuctuation of the particle 
number in cither pseudospin state of each species has a 
similar behavior as the entanglement, reaching the max- 
imal ^yN{N + 2)/12 as S" = [19], and reaching the 
minimum as S" = A^. The interspecies entanglement 
and polarization significantly affect the Gross-Pitaevskii 
equations governing the orbital wave functions associated 
with each pseudospin state of each species. This is a very 
interesting interplay between spin and orbital degrees of 
freedom. More phenomenology beyond that of the two- 
component BEG and can be experimentally observed is 
under investigation. 



Appendix A: and in all parameters regimes 
with B / 

In this and next appendices, we find S"" in all regimes 
of the parameters. Neighboring regimes with a same 
ground state can be combined, as described in the main 
text. We consider B in this Appendix. 
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FIG. 5. M2 = (Saz), for state \S, ±5), as a function of 5*. Each plot is for a given A'^ as indicated. 



1. B / 0, C > 



B ^0, C> 0, Ke = 



Were there not be the constraint ([T0|) on Sz , the mini- 
mum of the second term of ([S]) would be then at = Szo, 
where 

Szo = Int(-B/2C) 

is the integer closest to —B/2C. But one should, of 
course, consider the constraint pH)) . with the bounds ±5* 
determined by the sign of Ke- 



Now there is no Ke term. However, —Smax l£ S < 

I B I 

Smax is the largest possible range of S. If C < , 
we have S*™ — —sgn{B)Smax, which forces S"" — Smax- 
If <C < -J^, then SI" = Int(-2|) and thus 



Na+Nb — ^ — N^-N, 
|Int(-^)| < 5™ < ^^4^ 

B ^ and 



< S" 



If c > 

— 2 



\B\ 



then S" 



B^O,C>0,Ke>0 



a. B ^0, 00, Ke <0 

Then S^ = Smax, under which Sz can vary in the 
largest possible range. 5™ is determined by comparing 

SzO with the bounds ±Smax- If -Smax < SzO < Smax, 

that is, C > Tv^^W' ^^^"^ '^^^^ If 

iSzO ^ ~Smax, that is, C < jVa+A^i, ' tf^^'^ "^I" ~ —Smax- 

If 5*20 ^ Smax, that is, C < — j^J^j^^ , then S"™ = Smax- 



The minimum of the first term of ([81) is at 5 = Smin, 
which constrains the range of 5*2. As Sz term is positive, 
5™ may not be Smin- Nevertheless, if —Smin < Szo < 



-'mm, that is, C 2^ Na — Nb 
SzO- If •S'zO < —Smin Or 5*20 ^ "S", 
\B\ 



, then S"" = Smin and S*™ = 
mm I that is, if C < 



jr^zTTTb' then we have 5^" = -sgn(B)5'™. 

Sm will be determined below altogether for various pa- 
rameter regimes with S"™ = — sgn(_B)S'™. The present 
condition C > and K^ > overlaps only with the sub- 
case Ke > — C in Appendix lA 41 Hence there are three 
possibilities under the condition C < _ - 
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FIG. 6. The fluctuation ANca- , for state \S, ±5*), as a function of 5*. Each plot is for a given A'^ as indicated. 



2. C = 



In this case, in order that BSz is smallest, we have 
S™ = — sgn(i?)S'™, where sgn{B) represents the sign of 
B, S*™ is the value of S which minimizes E. Using the 
result of Appendix IA41 we know that there are three 
subcases. For Kg = 0, 5™ = Smax- For Ke > 0, the 
result of applies. For Ke < 0, the result of applies. 



3. B / 0, C < 

With C < 0, the S'^-dependent term in the energy 
dS]) is minimized always at 5*™ = — sgn(i?)S'm. This can 
be seen by representing this term as a parabola open- 
ing downward, with Szo the maximal point. All three 
subcases of Appendix apply under the condition C < 0. 



4. The cases with B / and & = -sgn(B)5 

As discussed above, for the case of B 7^ and C < 0, 
as well as the case of i? 7^ 0, C > while > 0, we 
have Sz = —sgn{B)S, therefore E can be expressed to 



E = {K, + C) S^ + {Ke-\B\)S 
^{Ke + C){S-Sof+E'i, 



where 



^0 



1^1 - 

2{Ke + C)'' 



(Al) 
(A2) 



(A3) 



is the the extreme point of the parabola {Ke + C){S 
S'o)^, Eq is independent of S and Sz- 



-C 



This condition have overlap with C < 0, while do not 
overlaps with the condition C > and Ke > 0. Then 



E— {Ke~\B\)S+const. Therefore we know the following 
subcases, which . 

liKe = -C < \B\ 

liKe = -C> \B\ 

liKe = -C= \B\ 

Q . < Q-m < C 



cm c 

cm c 

S*™ is any any integer in the range 



> ~C 
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2. B = 0, C = 



This condition overlaps with both the condition C < 
and the condition C > and Kf, > 0. 

Now that Ke + C >0, the parabola {Ke + C)(S'- S'o)^ 
open upwards. There are several subcases depending on 
the range of Sq, as one can conceive by considering the 
position of 5*0 related to the range (i) If ■S'o < Smin, 
which requires > ^^^N^^Nb+i'^ ' ^^^"^ ^"^ = Smin ■ (ii) 
If Smin < Sq < Smax, which rcquircs ^^^^^Jl'^^^^l'^ < 
Ke < which is a subset oi ~C < < 

\B\, then 5™ = So. (iii) If So > Smax, which requires 



C. Ke < -C 



There is no Sz term in E. Hence the values of S*™ 
are the same for C > 0, while 5™ can be any value 
in the legitimate ranges. If > 0, 5™ = Smin while 

-Sm^n < Sf < Sm^n■ K i^e < 0, 5™ = Smax while 
-Smax < Sf < Smax- liB = C = Ke = 0, thcu Sm and 

5™ are arbitrary values within the ranges ^ and (fTO|) . 

3. B = 0, C < 

Then we have Sz = iS in order to minimize E, con- 
sequently E can be written as Equations (jAip and (jA2l) , 
with _B = and thus 5*0 = — 2{k\c) • '^^^ discussions 
in App. lA 41 applv with B set to 0. But as B = causes 
both simplification and constraint. So we give details in 
the following. 



This condition overlaps with the condition C < 0, 
while do not overlap with the condition C > and 
Ke > 0. 

With Ke + C < 0, the coefficient in (|A2p is negative. 
The parabola {Kg + C){S — Sq)^ open downwards. 

We find that 5" = Smax for Ke < ^\+f , while 

S™ — Smin for TV +1 ^ ^C*, and S*™ = Smax 

and S*™ = 5'„ii„ are two degenerate solutions for Ke < 
■L^t^sii This result is obtained by considering the po- 
sition of 5*0 relative to the range ©, as the following, 
(i) So < Sa, which requires Ke < ^\+f , then S"" = 

Smax- (fi) So = Sa, which rcquircs Ke = '^j^+l'^ ' ^'^'^^'^ 
are two degenerate solutions S"" = S'maa; and S"" = Smin- 

(iii) When S'o > S'a, which requires Ke > ^^j^^l^ , then 
qni c ^ 

^ — '-^m.in. - 



a. B = 0, C < 0, = -C 

Now E = Xe-S* + const. If X,, -C < 0, S"" = Smax, 
__ If = (7 > 0, — Smin-; — 



b. B = 0, C <0, Ke> -C 

Now that Ke + C>0, the parabola {Ke + C){S - Sq)'^ 
open upwards. On the other hand now Ke > — C > 0, 



thus 5o = — 
-I- ■ 



< 0. Therefore S" 



q qm 



B ^0, C <0, Ke < -C 



Appendix B: S" 



and in all parameters regimes 
with 5 = 



Now we look at the special situation of S = 0. 



With Ke + C < 0, the parabola {Ke + C){S - So)'^ 
open downwards. We can directly use the result in 
App. IA4[ setting S = 0, to know that S"™ = Smax for 



Ke < 



NaC 



while S"- 



and S"" = Smax and 5"" 
solutions for Ke = — ^"^^ - 



S, 



min 

S; 



for 



Na + 1 < Ke < —C, 

are two degenerate 



1. B = 0, C > 

Then 5™ — in order to minimize CSz- Subsequently, 
it is easy to see that S*™ ~ Smin if Kg > 0, 5™ — Smax 
if Ke < 0, and 5™ can be any legitimate 5", i.e. Smin < 

gm ^ Smax , if Ke = . 
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